® REPUBLIQUE FRANQAISE 

INSTITUT NATIONAL 
DE LA PROPRIETE INDUSTRIELLE 



PARIS 



(0) N° de publication : 

(a n'utiliser que pour les 
commandes de reproduction) 

N° d enregistrement national 



(fl) IntCI 6 : H 04 M 9/08 



2 762 467 
97 04682 



® 



DEMANDE DE BREVET D'INVENTION 



A1 



(|2) Date de depot : 1 6.04.97. 
@) Priorite : 



Date de mise a la disposition du public de la 
demande : 23.10.98 Bulletin 98/43. 

@) Liste des documents cites dans le rapport de 
recherche pre Mm ina ire : Se reporter a la fin du 
present fascicule 

References a d'autres documents nationaux 
ap pa rentes : 



@)Demandeur(s) : FRANCE TELECOM SOCfETE ANO- 
NYM E — FR. 



@ Inventeur(s) : GILLOIRE ANDRE et TURBIN VALE- 
RIE. 



@)Titulaire(s) : 



© 



Mandataire(s) : CABINET PLASSERAUD. 



PROCEDE D'ANNULATION D'ECHO ACOUSTIQUE MULTI-VOIES ET ANNULEUR D'ECHO ACOUSTIQUE 
MULTI-VOIES. 

_ Pour effectuer une annulation d'echo sur plusieurs 
voies de signal sonore, sur chaque voie, on engendre un si- 
gnal synthetique ayant les caracteristiques spectrales d'un 
bruit blanc; on calcule des seuils de masquage frequentiel 
dans plusieurs bandes de frequence adjacentes, en utilisant 
des proprietes de la perception auditive; dans chaque ban- 
de de frequence, on amene le niveau du signal synthetique 
a la valeur du seuil de masquage associe, de facon a obtenir 
un signal auxiliaire, qu'on additionne au signal associe au 
haut-parleur de la voie consideree; on foumit le signal obte- 
nu a un filtre adaptatif; on additionne les signaux de sortie 
des filtres adaptatifs associes a chaque voie de micropho- 
ne; on soustrait le signal resultant du signal recu par le mi- 
crophone de cette voie; on en deduit une erreur 
d'estimation, en fonction de laquelle on corrige les coeffi- 
cients des filtres adaptatifs associes a cette voie. 





GENERATION DE 
SIGNAUX SYNTHETIQUES 










CALCUL DES 
SEUILS DE MASQUAGE 










GENERATION DE 
SIGNAUX AUXILIAIRES 








SUR CHAQUE VOIE, 
ADDITION DU SIGNAL AUXILIAIRE 
ET DU SIGNAL DE HAUT-PARLEUR 








FLLTRAGE ADAPTATIF 





PROPRIETES 
DE LA 

PERCEPTION 
AUDITIVE 



Hill 



2762467 



PROCEDE D'ANNULATION D 1 ECHO ACOUSTIQUE MULTI-VOIES ET 
ANNULEUR D T ECHO ACOUSTIQUE MULTI-VOIES 

La pr6sente invention concerne un proc6d£ d'annula- 
5 tion d'fecho acoustique multi-voies, ainsi qu'un annuleur 
d ' 6cho acoustique multi-voies . 

Dans le domaine de la transmission des signaux 
sonores, dans certaines applications telles que la t616pho- 
nie "mains libres" et la t616conf6rence, l ? 6cho acoustique 
10 est une source de g£ne considerable. Les dispositifs connus 
qui combattent l'fecho acoustique comportent gSneralement des 
f litres adaptatif s, dont le rdle est d' identifier et de 
mod61iser la reponse impulsionnelle du chemin de couplage 
acoustique entre le ou les haut-parleurs et le ou les 
15 microphones du syst&me de transmission de signaux sonores 
consid6r6 . 

La figure 1 rappelle la structure g6n6rale d f un 
annuleur d f 6cho acoustique classique. II est associ6 & un 
haut-parleur 10 et un microphone 12 entre lesquels il existe 

20 un chemin de couplage acoustique ou chemin d'6cho 14. Les 
r6f6rences 16 et 18 d&signent respectivement le signal 
sonore regu et le signal sonore emis. L f annuleur d'6cho de 
la figure 1 comprend un filtre adaptatif 20 recevant en 
entree le signal sonore regxx 16. L'objet du filtre adaptatif 

25 20 est d ! estimer par ses coefficients la r6ponse impulsion- 
nelle du chemin d'6cho 14, afin de soustraire l'echo du 
signal re<ju par le microphone 12. Pour cela, la sortie du 
filtre adaptatif 20 est connect6e & un soustracteur 22 qui 
soustrait le signal obtenu en sortie du filtre adaptatif 20 

30 du signal regu par le microphone 12. Le signal de difference 
obtenu en sortie du soustracteur 22 fournit une erreur 
d* estimation 24. Les coefficients du filtre adaptatif 20 
sont a Justus au cours du temps par un algorithme approprie 
qui exploite 1 ' information d ' erreur d 'estimation. 

35 Le choix de cet algorithme conditionne de fagon 

d6terminante les performances de 1' annuleur d'£cho. Les 
proc6d£s et dispositifs d'annulation d'Scho connus utilisent 
couramment un algorithme dit du gradient stochastique 
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normalise, habituellement d£sign6 par le sigle NLMS. Un 
inconvenient de cet algorithme est que sa vitesse de 
convergence depend des caractferistiques spectrales du signal 
sonore regu. Un objectif de la prSsente invention est de 
5 r&duire cette d6pendance. 

Le probl6me de 1 ' 6cho acoustique se pose bien sur 
egalement dans les systemes multi-voies de transmission de 
signaux sonores, c'est-d-dire dans les systemes d plusieurs 
haut-parleurs et plusieurs microphones, par exemple deux 
10 haut-parleurs et deux microphones dans le cas de la st6r6o- 
phonie. 

On connait des proc6d6s et des dispositifs d 1 ambula- 
tion d'6cho multi-voies fond6s sur le meme principe que ceux 
d'annulation mono-voie. La figure 2 repr6sente 4 titre 

15 d' exemple le schema de principe partiel d'un annuleur d'6cho 
acoustique stereophonique classique, ou, pour plus de 
clarte, on n'a represents qu'une seule des deux voies de 
microphone. La structure gfenferale peut fitre ais6ment 
g6n6ralis6e & celle d'un annuleur d'6cho acoustique £ N 

20 voies de signal sonore, ou N est un entier quelconque* On ne 
decrit dans ce qui suit que le cas N » 2. De fagon analogue 
a 1 1 annuleur d ' echo de la f igure 1 f 1 1 annuleur d 1 6cho 
st6r6ophonique s f applique & des voies de signal sonore regu 
16 lf 16 2 et des voies de signal sonore 6mis, dont une seule, 

25 d6sign6e par le chiffre de reference 18^ est representee* 
L' annuleur d'echo est associe & deux haut-parleurs 10^ 10 2 
et deux microphones, dont un seul, d6sign6 par le chiffre de 
reference 12 lf est represents. II existe quatre chemins 
d f 6cho, dont deux chemins d'6cho (14 x et 14 2 ) entre les deux 

30 haut-parleurs 10^ 10 2 et le premier microphone 12 : et deux 
autres chemins d'6cho (non repr6sent6s) entre les deux haut- 
parleurs 10 lf 10 2 et le second microphone. Pour estimer les 
reponses impulsionnelles des divers chemins d f 6cho, un 
filtre adaptatif est pr6vu entre chaque voie de haut-parleur 

35 et chaque voie de microphone. Ainsi, si on consid£re la voie 
du microphone 12 lf il est pr6vu deux filtres adaptatifs 20 lt 
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20 2 qui regoivent respect ivement en entr6e les signaux 
sonores regus et 16 2 . Les sorties des filtres adaptatifs 
20^ 20 2 sont fournies en entr6e A un additionneur 26 x . Dans 
un annuleur d'6cho k N voies ou N est un entier quelconque, 
5 on pr6voit un tel additionneur sur chaque voie de micropho- 
ne. La somme des signaux de sortie des filtres adaptatifs 
fournie par 1 1 additionneur 26 1 est soustraite, par un 
soustracteur 22 x , du signal regu par le microphone 12 x . La 
m£me op6ration est effectu6e sur chaque voie de microphone. 

10 Le signal de difference obtenu en sortie du soustracteur 22 L 
fournit une erreur d' estimation 24 x commune d tous les 
filtres adaptatifs relatifs a la voie de microphone considfe- 
r6e ; dans l'exemple de la figure 2, l f erreur d' estimation 
24 x est commune aux filtres adaptatifs 20 x et 20 2 . De mgme 

15 qu'en annulation d f 6cho mono- voies, les coefficients des 
filtres adaptatifs sont modifies de fagon iterative par un 
algorithme appropri6, sur chaque voie de microphone, & 
partir de 1' erreur d ? estimation obtenue. 

On a constats, en annulation d'6cho multi-voies, et 

20 en particulier en annulation d'6cho ster6ophonique, que la 
vitesse de convergence des algorithmes d* adaptation des 
coefficients des filtres adaptatifs etait comparat ivement 
plus faible que la vitesse de convergence de ces algorithmes 
dans le cas de l f annulation d'6cho mono-voie. II a 6t6 

25 montr§ que ce ralentissement de la convergence provient du 
fait que les signaux sonores regus par les haut-parleurs, 
d£sign£s par les chiffres de reference 16 x et 16 2 sur la 
figure 2, sont mutuellement corr616s. 

Le ralentissement de la convergence entraine 

30 plusieurs inconvfenients. En particulier, dans un systdme de 
t616conf6rence, les locuteurs situfes dans la salle distante 
pergoivent d 1 autant plus longtemps 1 f 6cho de leur parole 
lors de la mise en route du syst&me ou aprSs une modifica- 
tion acoustique (mouvements des auditeurs par exemple) dans 

35 la salle ou est situe 1 1 annuleur d'6cho. En outre, en 
annulation d'6cho multi-voies, toute modification acoustique 
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dans la salle distante perturbe la convergence des filtres 
adaptatifs de l'annuleur d f 6cho, en raison de la correlation 
mutuelle, 6voqu6e pr6c6demment , entre les signaux de parole 
regue, ce qui entraine une r6apparition ou une augmentation 
5 du niveau de 1 1 echo . 

En revanche, on a observe en pratique que lorsque 
des composantes non mutuellement corr616es sont pr6sentes 
sur chaque signal microphonique, el les tendent & acc616rer 
la convergence des annuleurs d*6cho multi-voies* Un objectif 

10 de la prSsente invention est d 1 exploiter cette propri6t6 des 
composantes mutuellement d6corr616es pour am61iorer les 
performances des annuleurs d'6cho multi-voies, et am&liorer 
par consequent la qua lite de la communication dans les 
systSmes de transmission de signaux sonores qui mettent en 

15 oeuvre des annuleurs d ! 6cho multi-voies. Pour cela, le 
principe g6n6ral de la prfesente invention consiste 4 a j outer 
aux voies de signal sonore re$u des signaux auxiliaires 
mutuellement decorr616s, rendus inaudibles par exploitation 
de certaines propri6t6s de 1' audition humaine. 

20 Plus pr6cisement, la pr6sente invention propose un 

proc£d£ d'annulation d f 6cho sur N voies de signal sonore 
ayant chacune un haut-parleur et un microphone associ6s, N 
6tant un entier sup6rieur ou 6gal k 1, suivant lequel, sur 
chacune des N voies : 

25 (a) on engendre un signal synthetique ayant les 

caracteristiques spectrales d ! un bruit blanc, le spectre de 
ce signal s'etendant sur plusieurs bandes de frequence 
adjacentes, et ce signal synthetique fetant d6corr616 des 
signaux synthetiques engendr6s sur les autres voies ; 

30 (b) pour chaque bande de frequence, on calcule un 

seuil de masquage frfequentiel correspondant au signal 
associS au haut-parleur de la voie consid6r6e en utilisant 
des propri6t6s de la perception auditive humaine ; 

(c) dans chaque bande de frequence, on am&ne le 

35 niveau du signal synthetique d la valeur du seuil de 
masquage fr6quentiel associe, de fagon & obtenir un signal 
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auxiliaire ; 

(d) on additionne le signal auxiliaire au signal 
associ6 au haut-parleur de la voie considferfee, les seuils de 
masquage fr6quentiel ayant 6t6 pr6c6demment calculus de 

5 fagon d rendre le signal auxiliaire inaudible, et les 
signaux auxiliaires des N voies 6tant mutuellement d6corr6- 
16s ; 

(e) on fournit le signal obtenu pr6c6demment en 
entree & un filtre adaptatif dont les coefficients forment 

10 une estimation de la r£ponse impulsionnelle du chemin de 
couplage acoustique entre le haut-parleur et le microphone 
associ&s d la voie de signal sonore consid£r£e ; 

( f ) on additionne les signaux obtenus en sortie des 
N filtres adaptatifs respectivement associ&s & chaque voie 

15 de microphone et on soustrait le signal resultant du signal 
regu par le microphone assocife & cette voie ; 

(g) on calcule une erreur d 1 estimation d partir de 
la difference obtenue h 1' issue de la soustraction pr6c6- 
dente ; 

20 (h) on corrige les coefficients des filtres adapta- 

tifs associfes a la voie de microphone consid6r6e en fonction 
de l r erreur d f estimation associ&e. 

Dans un mode particulier de realisation du proc6d£, 
pour calculer chaque seuil de masquage fr6quentiel, 
25 (bl) on forme des blocs contenant chacun un nombre 

pr6d6termin6 d'6chantillons de signal sonore, deux blocs 
successifs se recouvrant mutuellement sur un nombre prede- 
termine d 1 fechanti lions ; 

(b2) on pond&re les 6chantillons de chaque bloc par 
30 une fenfitre d ! apodisation ; 

(b3) on calcule la transform6e de Fourier de chaque 

bloc ; 

(b4) on divise le domaine des frequences en plu- 
sieurs bandes critiques adjacentes ayant des frequences 
35 infferieures et sup&rieures determines ; 

puis, pour chaque bloc et dans chaque bande critique : 
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(b5) on additionne l'6nergie des differentes raies 
spectrales appartenant d la bande critique, de fagon d 
obtenir la valeur, pour cette bande critique, d'une fonction 
de repartition de l'6nergie ; 
5 (b6) on calcule le produit de convolution de la 

fonction de repartition de l'6nergie et d'une fonction 
d'6talement basilaire, obtenue d partir d'une table de 
consultation, de fagon d obtenir un spectre 6tal6 E t ; 

(b7) on calcule un indice de tonalite a L d partir des 
10 energies des raies spectrales appartenant aux differentes 
bandes critiques ; 

(b8) on calcule un facteur de correction 0 4 & partir 
de 1' indice de tonalite calcule pr6c6demment ; 

(b9) on calcule un seuil de masquage fr6quentiel T t 
15 & partir du spectre etaie E L et du facteur de correction 0 if 
d6fini comme suit : 

10 log 10 T 4 = 10 log 10 E t - 0 t - C t 
ou C 1 est un parametre suppl6mentaire de correction de seuil 
calculi S partir de plusieurs caract6ristiques des signaux 
20 sonores des N voies. 

Dans un mode particulier de realisation du proc6d6, 
les caract6ristiques servant d calculer le parametre 
suppl6mentaire de correction de seuil C t comprennent le 
niveau respect if des signaux des N voies et un ensemble de 
25 valeurs de la fonction d ' intercorr61ation d'au moins 
certaines paires de ces signaux, le parametre suppl6mentaire 
de correction C L 6tant une fonction croissante des valeurs 
de cette fonction d'intercorr61ation. 

En variante, les caract6ristiques servant £ calculer 
30 le parametre supplemental re de correction de seuil C : 
peuvent comprendre la somme des Energies contenues dans les 
diver ses bandes critiques de chaque voie et le rapport des 
Energies par bande critique d'au moins certaines paires des 
differentes voies, le parametre supplement aire de correction 
35 C t etant d'autant plus important que la somme des energies 
contenues dans les diverses bandes critiques est import ante 
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et que le rapport des energies par bande critique d'une 
paire de voles est peu different entre les diff£rentes 
bandes critiques. 

La pr6sente invention propose 6galement un annuleur 
5 d f 6cho & N voies de signal sonore ayant chacune un haut- 
parleur et un microphone associfes, N 6tant un entier 
supferieur ou 6gal & 1, comprenant : 

- des moyens de calcul de signaux auxiliaires, 
comport ant : 

10 des moyens de formation de blocs contenant 

chacun un nombre pr6d6termin6 d'6chanti lions du 
signal sonore, deux blocs successifs se recouvrant 
mutuellement sur un nombre predetermine d ! 6chantil- 
lons ; 

15 des moyens de pond6ration des 6chantillons de 

chaque bloc par une fen£tre d f apodisation, disposes 
en sortie des moyens de formation de blocs ; 

des moyens de calcul de transformee de 
Fourier, disposes en sortie des moyens de pond6ra- 

20 tion ; 

des moyens de calcul des valeurs de 
frequences infferieures et sup6rieures de bandes 
critiques adjacentes, ou une table de consultation 
contenant ces valeurs ; 

25 des moyens pour calculer une fonction de 

repartition de l'6nergie sur les bandes critiques, 
la valeur, dans une bande critique donn6e, de la 
fonction de repartition de l'6nergie 6tant calcul6e 
en additionnant l'6nergie des diff6rentes raies 

30 spectrales appartenant a cette bande critique ; 

une seconde table de consultation, contenant 
des valeurs repr6sentant des fonctions d f 6talement 
basilaire associ&es chacune & une bande critique 
donn6e ; 

35 des moyens pour calculer un spectre etaie E 1# 

en calculant, pour une bande critique donn6e, le 
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produit de convolution de la fonction de r6partition 
de l'&nergie et de la fonction d'&talement assocides 
& cette bande critique ; 

des moyens de calcul d'un facteur de correc- 
5 tion 0 t associ& & une bande critique donn6e ; 

des moyens de calcul d'un seuil de masquage 
frfequentiel T L associ6 k une bande critique donn6e, 
d partir du spectre 6tal6 et du facteur de correc- 
tion 0 i associ6s £ cette bande critique, le seuil de 
10 masquage frfequentiel T t 6tant defini comme suit : 

10 log 10 Ti = 10 log 10 E L - 0 ± - Ci 
ou C L est un paramdtre suppl&nentaire de correction 
de seuil ; 

des moyens de g6n6ration de signaux synth6ti- 
15 ques ayant les caract6ristiques spectrales d'un 

bruit blanc ; 

des moyens de contrdle, connect6s aux N voies 
de signal sonore, pour calculer pour chaque bande 
critique le param6tre supplfementaire de correction 
20 de seuil C x S partir de plusieurs caract6ristiques 

des signaux sonores des N voies ; 

des moyens pour amener le niveau d'un signal 
synthfetique a la valeur du seuil de masquage frfe- 
quentiel associ6, de fa<jon a obtenir un signal 
25 auxiliaire ; 

des moyens de calcul de transform6e de Fourier 
inverse ; 

des moyens de mfemoire pour stocker le r6sultat 
du traitement des blocs successifs ; 
30 des moyens pour combiner le rdsultat du 

traitement d'un bloc avec le r6sultat du traitement 
du bloc pr6c6dent ; 

des moyens de lecture sequentielle reli&s 
auxdits moyens de m&noire, pour faire passer la 
35 cadence des donnSes de la cadence bloc & la cadence 

§chantillon ; 
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N premiers moyens additionneurs respectivement 
disposes sur les N voies de signal sonore, pour 
additionner sur chaque voie le signal sonore re<ju et 
le signal auxiliaire associ6 ; 
5 - NxN moyens de filtrage adaptatif, comprenant N 

moyens de filtrage adaptatif associ&s & chacune des N voies 
de microphone, et ayant chacun des coefficients qui forment 
une estimation de la r£ponse impulsionnelle du chemin de 
couplage acoustique entre le haut-parleur et le microphone 
10 associ6s & une des N voies de signal sonore ; 

- des seconds moyens additionneurs, connect6s en 
sortie des N moyens de filtrage adaptatif associds d chaque 
voie de microphone, pour additionner les signaux de sortie 
de ces N moyens de filtrage adaptatif ; 

15 - N moyens soustracteur s , disposes chacun en sortie 

des seconds moyens additionneurs et connectds respectivement 
aux microphones des N voies de signal sonore, pour sous- 
traire sur chaque voie le signal obtenu en sortie des 
seconds moyens additionneurs du signal regu par le micro- 

20 phone de cette voie ; 

- N moyens pour calculer sur chaque voie une erreur 
d' estimation £ partir du r6sultat fourni par le moyen 
soustracteur connects au microphone de cette voie ; 

- N moyens pour corriger de fa<jon iterative les 
25 coefficients respectifs des N moyens de filtrage adaptatif 

associ&s & chacune des N voies de microphone, en fonction de 
1 ' erreur d f estimation associfee & chaque voie de microphone. 

Les realisations different selon les applications, 
et en particulier, le nombre de haut-parleurs du syst6me de 

30 transmission de signaux sonores consider^ peut etre diffe- 
rent du nombre de microphones • 

Dans un mode particulier de realisation du disposi- 
tif, les caract6ristiques servant S calculer le paramdtre 
suppl6mentaire de correction de seuil Ci comprennent le 

35 niveau respectif des signaux des N voies et un ensemble de 
valeurs de la fonction d' intercorrelation d'au moins 
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certaines paires de ces signaux, le param6tre suppl&mentaire 
de correction C A 6tant une fonction croissante des valeurs 
de cette fonction d' intercorr61ation. 

En variante, les caract6ristiques servant & calculer 
5 le paramdtre supplSmentaire de correction de seuil C t 
peuvent comprendre la somme des Energies contenues dans les 
diverses bandes critiques de chaque voie et le rapport des 
fenergies par bande critique d'au moins certaines paires des 
diffSrentes voies, le param&tre suppl6mentaire de correction 

10 C i 6tant d ' autant plus important que la somme des Energies 
contenues dans les diverses bandes critiques est importante 
et que le rapport des Energies par bande critique d v une 
paire de voies est peu different entre les diff£rentes 
bandes critiques. 

15 Dans un mode particulier de realisation du disposi- 

tif, l'annuleur d'6cho comprend en outre des moyens de 
calcul d'un indice de tonalit6 a t & partir des Energies des 
raies spectrales appar tenant aux diffSrentes bandes criti- 
ques, lesdits moyens de calcul du facteur de correction 0 t 

20 calculant le facteur de correction 0 ± & partir de 1 1 indice 
de tonalit6 a A . 

Dans un mode particulier de r6alisation du disposi- 
tif, les moyens de calcul de 1 1 indice de tonalit6 fournis- 
sent un indice de tonalite constant, identique pour toutes 

25 les bandes critiques et d6fini par : 

a, = a = min(SFM/SFM aax/ 1) 
ou min(a, b) d&signe la plus petite des valeurs a et b, 
ou SFM^ est un paramdtre de valeur pr&d&termin&e en dB 
associS 4 un signal sinusoidal pur, et 

30 ou SFM « 10 log 10 G/A, ou log 10 d6signe le logarithme dfecimal, 
G d6signe la moyenne g6om6trique de l'6nergie sur un nombre 
pr6d6 termini de points de la trans formee de Fourier et A 
designe sa moyenne arithm6tique sur le mdme nombre de 
points - 

35 

Dans un mode particulier de rfealisation, SFM^ * 
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- 60 dB. 

Dans un mode particulier de realisation, les moyens 
de calcul du facteur de correction 0 t fournissent un facteur 
de correction 0 L d6fini par : 
5 0i - max (SO, a L . (k x + B t ) + (1 - a^.kj) 

ou max(a, b) d6signe la plus grande des valeurs a et b, 

ou SO, k x et k 2 sont des param6tres de valeur pr6d6termin6e 

en dB, 

ou a t est l'indice de tonalite associe & la bande critique 

10 consideree, et 

ou B t d6signe la frequence de la bande critique en Bark. 

Dans un mode particulier de realisation du disposi- 
tif, SO » 24,5 dB, k x = 14,5 dB et k 2 - 5,5 dB. 

D'autres caract6ristiques et avantages de la 

15 pr&sente invention apparaitront & la lecture de la descrip- 
tion d6taill6e qui suit de modes particuliers de realisa- 
tion, donnfes A titre d'exemples non limitatifs. La descrip- 
tion se r6f6re aux dessins annexes, dans lesquels : 

- la figure 1, relative & I'art ant6rieur et d6crite 
20 en introduction, est une representation schematique d ? un 

annuleur d f 6cho acoustique mono-voie classique ; 

- la figure 2, 6galement relative & l'art anterieur 
et d6crite en introduction, est une representation schemati- 
que partielle d'un annuleur d'6cho acoustique st6r6ophonique 

25 classique, ou une seule voie de microphone est representee ; 

- la figure 3 est un organigramme d6crivant la 
succession g6n6rale des etapes du proc6d6 de la pr6sente 
invention dans un mode particulier de realisation ; 

- la figure 4 repr6sente de fagon sch6matique le 
30 dispositif de la pr6sente invention dans un mode particulier 

de realisation ; 

- la figure 5 repr6sente de fagon sch6matique la 
structure du module de calcul de signaux auxiliaires compris 
dans le dispositif de la pr6sente invention, dans un mode 

35 particulier de realisation ; 
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- la figure 6a est un graphe illustrant un premier 
exemple devolution au cours du temps de l'fecart relatif 
entre les r&ponses impulsionnelles r^elles des chemins de 
couplage acoustigue et les r6ponses impulsionnelles estimates 

5 par les f litres adaptatifs d'un annuleur d'Scho st6r6ophoni- 
que conforme £ un mode particulier de realisation de la 
pr6sente invention ; 

- la figure 6b est un graphe analogue & celui de la 
figure 6a, etabli pour un annuleur d'6cho st6r6ophonique 

10 classique ; 

- la figure 7a est un graphe illustrant un deuxidme 
exemple devolution au cours du temps de l'6cart relatif 
entre les r6ponses impulsionnelles r belles des chemins de 
couplage acoustique et les r6ponses impulsionnelles estimdes 

15 par les f litres adaptatifs d'un annuleur d'6cho st6r6ophoni- 
que conforme & un mode particulier de realisation de la 
prfesente invention ; 

- la figure 7b est un graphe analogue d celui de la 
figure 7a, 6tabli pour un annuleur d f 6cho st6r6ophonique 

20 classique* 

Le proc&d& et le dispositif de la prfesente invention 
s'appliquent & un signal sonore quelconque. A titre d' exem- 
ple non limitatif, il peut s'agir aussi bien d'un signal de 
parole que de musique. 

25 On va d'abord d6crire 1 'organisation g£n£rale du 

proc£d£ d'annulation d'£cho multi-voies dans un mode 
particulier de realisation, & 1 9 aide de la figure 3. 

On considdre N voies de signal sonore dont chacune 
comporte un haut-parleur et un microphone. 

30 Sur chaque voie de signal sonore regu, on synth6tise 

tout d'abord un signal ayant les caract£ristiques spectrales 
d'un bruit blanc et dont le spectre s f 6tend sur plusieurs 
bandes de frequence adjacentes, les signaux synthfetiques des 
N voies 6tant mutuellement decorr616s. 



2762467 



13 

Ensuite, pour chaque bande de frequence, on calcule 
un seuil de masquage frequentiel correspondant au signal 
sonore regu par le haut-parleur de la voie consid6r6e. 

On rappelle que l'effet de masquage consiste en ce 
5 qu'un son en presence d'un autre peut devenir partiellement 
ou compietement inaudible. Lorsque ce ph6nom6ne se produit 
dans le domaine fr6quentiel, on parle de masquage simultanfe 
ou masquage fr6quentiel. II existe un seuil, appeie seuil de 
masquage fr6quentiel, en de<?& duquel toutes les composantes 

10 frequentielles d f un son masqu6 sont inaudibles. Pour plus de 
detail sur le ph6nom6ne de masquage, on se reportera 
utilement £ l'ouvrage de E * ZWICKER et R. FELDTKELLER 
intitule "Das Ohr als Nachrichtenemp f&nger " , Stuttgart, West 
Germany, Hirzel Verlag, 1967. 

15 Dans un mode particulier de realisation, le calcul 

des seuils de masquage frequent iel s'effectue comme suit. 
Tout d'abord, on transforme la cadence du flot de donn6es de 
signal sonore en passant d'une cadence d'6chanti lions, par 
exemple de 16 kHz, a une cadence de blocs, par exemple 

20 chacun d ? une dur6e de 32 ms. Chaque bloc contient done un 
nombre predetermine d f 6chantillons . Deux blocs successifs se 
recouvrent mutuellement sur un nombre predetermine d f 6chan- 
tillons, par exemple sur une duree de 16 ms. 

Ensuite, pour am61iorer l'estimation ulterieure du 

25 spectre, on pond6re les echanti lions de chaque bloc par une 
fenStre d f apodisation, par exemple du type fenfitre de 
Hamming ou de Hanning. 

Puis on calcule la transform6e de Fourier de chaque 
bloc sur un nombre predetermine de points, par exemple 512 

30 points pour des blocs de 32 ms. On peut utiliser a cet effet 
un algorithme de transform6e de Fourier rapide. 

On divise ensuite le domaine des frequences auquel 
le signal sonore regu est susceptible d'appartenir en 
plusieurs bandes de frequences adjacentes ayant des fr6quen- 

35 ces inf6rieures et superieures d6termin6es. Ces bandes de 
frequences sont appeiees bandes critiques. Le tableau ci- 
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dessous fournit, & titre d'exemple non limitatif, la 
decomposition du domaine des frequences de la parole en 
bandes critiques (bornes infferieures et supferieures et 
frequences centrales) etablie par E . ZWICKER et R. FELDTKEL- 
5 LER et figurant dans l'ouvrage pr6cit6. La largeur de chaque 
bande critique est par convention egale & 1 Bark. Le tableau 
donne done la correspondance entre l'6chelle f r6quentielle 
des Hertz et celle des Bark. 

10 



40 



Dan /4 a 

critique 
n° 


inf&rieure 
(Hz) 


Centre 
(Ha) 


Borne 
sup&rieure 
(Hz) 




o 


50 


100 


2 


100 


150 


200 


3 


200 


250 


300 


4 


300 


350 


400 


5 


400 


450 


510 


6 


510 


570 


630 


7 


630 


700 


770 


8 


770 


840 


920 


9 


920 


1000 


1080 


10 


1080 


1170 


1270 


11 


1270 


1370 


1480 


12 


1480 


1600 


1720 


13 


1720 


1850 


2000 


14 


2000 


2150 


2320 


15 


2320 


2500 


2700 


16 


2700 


2900 


3150 


17 


3150 


3400 


3700 


18 


3700 


4000 


4400 


19 


4400 


4800 


5300 


20 


5300 


5800 


6400 


21 


6400 


7000 


7700 


22 


7700 


8500 


9500 


23 


9500 


10500 


12000 


24 


12000 


13500 


15500 



Au lieu d'appliquer une transform6e de Fourier conune 
d6crit ci-dessus, on peut, en variante, utiliser une 
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decomposition par bancs de f litres, ce qui permet d'obtenir 
directement la decomposition en bandes critiques. 

L'6tape suivante consiste & additionner l f 6nergie 
des dif ffirentes raies spectrales appartenant d chaque bande 
5 critique, de fagon a obtenir la valeur, dans cette bande 
critique, d'une fonction de repartition de l'6nergie. 

On calcule ensuite le spectre etaie E if dfefini comme 
le produit de convolution de la fonction de repartition de 
l f 6nergie par bande critique avec la fonction dite d'6tale- 

10 ment basilaire. La fonction d'6talement basilaire E b 
correspond S la r6ponse f r6quentielle en 6nergie de la 
membrane basilaire de l'oreille humaine. Elle repr6sente 
l f excitation resultant de la perception d'une tonale, c'est- 
d-dire d'un son de frequence pure, et peut etre definie par 

15 1' equation suivante, figurant 6galement dans l'ouvrage 
pr6cit6 de E. ZWICKER et R. FELDTKELLER : 
10.1og 10 E b = 15,81+7, 5. (B+0,474)-17, 5. ( l + (B+0,474) 2 ) l/2 
oH log 10 repr6sente le logarithme decimal et B est la 
frequence en Bark. Dans le mode particulier de realisation 

20 d6crit, les valeurs de cette fonction sont lues dans une 
table de consultation. 

On a montre experimentalement qu'il est plus 
difficile de masquer un bruit par une tonale que 1* inverse. 
Pour tenir compte de cette asym6trie des seuils de masquage, 

25 on applique k chaque seuil un facteur de correction 0 L 
associ6 & une bande critique donn6e B £ et qui depend de 
1' aspect tonal ou non tonal du signal sonore. Pour 6valuer 
cet aspect tonal ou non tonal, on utilise, dans un mode 
particulier de realisation, une mesure notee SFM (en anglais 

30 Spectral Flatness Measure), d6finie comme le rapport entre 
la moyenne g6om6trique de l'energie sur un nombre predeter- 
mine de points de la transform6e de Fourier et sa moyenne 
arithmetique sur le meme nombre de points, par exemple 256 
points pour une transform6e de Fourier calcuiee sur 512 

35 points. 
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La valeur SFM obtenue est compar&e k celle d'un 
signal sinusoidal pur, pour lequel, & titre d'exemple non 
limitatif , on peut poser par convention SFM • SFM^ = 
- 60 dB. 

5 On calcule alors un indice de tonalit6 a if qui, dans 

un mode particulier de realisation, est constant et identi- 
que pour toutes les bandes critiques et dfefini par : 

a t = a = min (SFM/SFM^,, 1) 
ou min (a, b) dfesigne la plus petite des valeurs a et b. 
10 L' indice a A est compris entre 0 et 1 et est plus 61ev6 pour 
les composantes tonales du signal que pour les composantes 
non tonales. Toute autre mfethode de calcul de 1' indice de 
tonalitfe jugee approprifee est envisageable. 

Le facteur de correction 0 t assocife A la bande 
15 critique de frequence B t en Bark est d6fini a partir de 
1' indice de tonalit6 a L comme suit : 

O i - max (SO, a 1 .(k l + BJ + (1 - a L ) .k 2 ) 
ou max (a, b) dfesigne la plus grande des valeurs a et b et 
ou SO, k x et k 2 sont des paramfetres de valeur pr6d6termin6e 
20 en dB. Typiquement, on peut choisir SO « 24,5 dB, k x « 14,5 
dB et k 2 « 5, 5 dB. 

Le facteur de correction 0 L est une fonction crois- 
sante du rang de la bande critique dans le cas de composan- 
tes tonales dominantes et accorde plus de poids & celles-ci 
25 qu'aux composantes non tonales. 

On applique en outre aux seuils de masquage un 
paramStre supplfementaire de correction C A calculfe & partir 
de plusieurs caract6ristiques des signaux sonores des N 
voles. Le paramfetre C A peut etre commun & tous les seuils de 
30 masquage correspondant aux diffSrentes bandes critiques, ou 
bien avoir une valeur sp6cifique pour chaque bande critique. 

Dans un mode particulier de realisation, on calcule 
le param^tre supplfementaire de correction de seuil C t & 
partir des niveaux des signaux sonores regus sur les N voies 
35 et d'un ensemble de valeurs de la fonction d' intercorrfela- 
tion d'au moins certaines paires de ces signaux, Dans ce 



2762467 



17 

mode particulier de realisation, le parametre suppl6mentaire 
de correction de seuil C A est une fonction croissante des 
valeurs de la fonction d' intercorr61ation. 

A titre d'exemple non limitatif, dans une installa- 
5 tion pratique typique de t616conf6rence h plusieurs micro- 
phones et plusieurs haut-parleurs, les microphones sont 
disposes de telle sorte que la correlation entre signaux 
captfes est importante entre deux microphones adjacents et 
faible entre microphones distants. Cette correlation 

10 variable se retrouve, moyennant des modifications 6ventuel- 
les dues aux traitements durant la transmission, d la 
reception, entre les voies adjacentes d'une part, & savoir, 
entre deux haut-parleurs directement voisins, et entre les 
voies distantes d' autre part, & savoir, entre les autres 

15 paires de haut-parleurs distants. 

Ainsi, on peut se limiter & utiliser, pour le calcul 
des paramdtres C if les caract6ristiques de niveau et 
d f intercorr61ation des paires de voies correspondant aux 
paires de haut-parleurs directement voisins et aux deux 

20 paires adjacentes (si ces derniferes existent, compte tenu du 
nombre total de haut-parleurs). Ainsi, si on considfere une 
paire de haut-parleurs directement voisins (HP k , HP k+1 ), on 
utilisera leurs caracteristiques de niveau et d ' intercorr6- 
lation ainsi que celles des paires adjacentes (HP k _ x , HP k ) 

25 (HPj,.! et HP k 6tant directement voisins) et (HP k+1 , HP k+2 ) (HP k+1 
et HP k+2 6tant directement voisins). 

Dans un autre mode de realisation, on calcule le 
parametre supplemental re de correction C L & partir de la 
sonune des energies contenues dans les diverses bandes 

30 critiques de chaque voie et le rapport des energies par 
bande critique d f au moins certaines paires des diff6rentes 
voies • 

De fa<jon analogue au mode de realisation pr6c6dent, 
on peut se limiter £ utiliser les caracteristiques de sonune 
35 et de rapport d'energie des paires de voies correspondant 
aux paires de haut-parleurs directement voisins et aux deux 
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paires adjacentes. Le paramdtre C t est d'autant plus 
important que la sonune des energies contenues dans les 
diverses bandes critiques est importante et que le rapport 
des Energies par bande critique d'une paire de voies est peu 
5 different entre les diff6rentes bandes critiques. 

Ainsi, lorsque 1 1 intercorrfelation entre les signaux 
sonores de deux voies donn6es est forte, et/ou que ces deux 
signaux ont des energies comparables, le paramfrtre C £ est 
61ev6, c'est-^-dire qu'une reduction importante du seuil est 

10 appliqu6e, typiquement de 6 4 10 dB. A 1' inverse, lorsque 
les signaux sont peu correl6s et/ou sont d' energies sensi- 
blement diff6rentes, une reduction moins importante du seuil 
est appliqu6e, typiquement de 0 & 6 dB, Cette correction 
tient compte des propri6t6s de 1* audition. Elle a pour but 

15 de rendre inaudibles les signaux auxiliaires mentionn6s en 
introduction et dont l'obtention est d6taill6e ci-aprSs, 
tout en permettant de leur donner le niveau le plus 61ev6 
possible. 

Pour chaque bloc et dans chaque bande critique, on 

20 obtient le seuil de masquage fr&quentiel T t & partir du 
spectre 6tal6 E t et du facteur de correction 0i definis 
precedemment, et du paramfetre supplementaire de correction 
de seuil C if comme suit : 

10 log 10 T L = 10 log 10 E L - 0 t - C t 

25 Comme le montre la figure 3, l'6tape suivante du 

procede consiste h obtenir, sur chaque voie, un signal 
auxiliaire. Pour cela, dans un mode particulier de realisa- 
tion, on amfene dans chaque bande critique le niveau du 
signal synthetique 5 la valeur du seuil de masquage frfequen- 

30 tiel associ6. 

On calcule ensuite la transform6e de Fourier inverse 
correspondant au r6sultat du traitement du bloc courant, par 
exemple au moyen d'un algorithme de transformfee de Fourier 
inverse rapide. Dans le cas ou on n'a pas utilise une 

35 transform6e de Fourier, mais une decomposition par bancs de 
f litres, comme d6crit pr6c6demment , on effectue f au lieu 



2762467 



19 

d'une transformation de Fourier inverse, une &tape de 
synthase correspondent S une transformation inverse de la 
decomposition par bancs de f litres. 

On combine ensuite le resultat du traitement des 
5 blocs precedents avec le r6sultat du traitement du bloc 
courant, afin de reconstruire le bloc de bruit masque 
courant, par exemple au moyen d'un proc6d6 connu de synthase 
par recouvrement et addition (en anglais " overlap- add H ) ou 
par recouvrement et sauvegarde (en anglais " overlap- save" ) . 
10 On repasse ensuite de la cadence de blocs d la 

cadence d'echantillons, 

L f ensemble des traitements d6crits pr6c6demment peut 
etre effectu6 par un processeur de traitement de signal en 
temps r6el, 

15 Dans un mode particulier de realisation, l'etape 

suivante du proc&d6 consiste, sur chaque voie, & additionner 
le signal auxiliaire obtenu au signal associe au haut- 
parleur de cette voie, comme l'indique la figure 3. Par 
construction, les signaux auxiliaires des N voies sont 

20 mutuellement d6corr616s et, du fait que les signaux auxi- 
liaires ne d6passent pas les seuils de masquage fr6quentiel 
associ6s, ils sont inaudibles. 

Dans le mode de realisation d&crit, chaque signal 
auxiliaire est additionn6 au signal de haut-parleur associe 

25 apr6s avoir 6t6 reconstitu6 dans sa totalite & partir des 
differentes bandes critiques. En variante, cette addition 
peut tout aussi bien Stre effectu6e dans chacune des 
diff6rentes bandes critiques. 

On effectue ensuite une etape globale de filtrage 

30 adaptatif, comme l'indique la figure 3. Cette 6tape est 
d6taill6e ci-dessous dans un mode particulier de realisa- 
tion • 

On fournit tout d'abord la somme du signal auxi- 
liaire et du signal de haut-parleur en entr6e & un filtre 
35 adaptatif dont les coefficients forment une estimation de la 
r&ponse impulsionnelle du chemin de couplage acoustique 
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entre le haut-parleur et le microphone associ6s & la voie de 
signal sonore consid6r6e. 

On procfede ainsi sur chaque voie, puis on additionne 
les signaux obtenus en sortie des N f litres adaptatifs 
5 associ6s a chaque voie de microphone. On soustrait ensuite 
le signal resultant, sur chaque voie de microphone, du 
signal req;u par le microphone associ6 a cette voie. 

A partir du signal de difference ainsi obtenu, on 
calcule sur chaque voie de microphone une erreur d'estima- 
10 tion et on corrige les coefficients des filtres adaptatifs 
associ&s A cette voie de microphone en fonction de cette 
erreur d' estimation. 

L'6tape de filtrage adaptatif d6taill6e ci-dessus 
est r6alis6e dans le domaine temporel. En variante, on peut 
15 rfealiser le filtrage adaptatif par blocs et non par 6chan- 
tillons, en particulier en effectuant 1' ensemble des traite- 
ments correspondants dans le domaine de Fourier. Pour une 
description de proc6d6s de filtrage adaptatif fr6quentiel, 
on se reportera utilement & 1' article de J. PRADO et E. 
20 MOULINES intitule "Frequency domain adaptative filtering 
with applications to acoustic echo cancellation", Annales 
des Telecommunications , 49, n° 7-8, 1994, pages 414-428. 

On d6crit dans ce qui suit un mode particulier de 
realisation d'un annuleur d'6cho conforme h la pr6sente 
25 invention, A l'aide des figures 4 et 5. 

Dans ce mode particulier de realisation, le nombre 
de haut-parleurs est identique au nombre de microphones. 
Cependant, diverses realisations sont possibles selon les 
applications, et en particulier, le nombre de haut-parleurs 
30 peut Stre different du nombre de microphones. 

Comme le montre la figure 4, 1' annuleur d'6cho est 
associe £ N haut-parleurs 30^ 30 H et N microphones 32 x , 

32 K . II comprend un module 34 de calcul de signaux 
auxiliaires dont la structure detaill6e sera d6crite & 
35 l'aide de la figure 5. De fa$on globale, le module 34 de 
calcul de signaux auxiliaires comprend un module de contrOle 
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36 connecte aux N voies de haut-parleur , N modules 38 lf . .., 
38 N de production de slgnaux auxiliaires respectivement 
connectfes aux N voies de haut-parleur et au module de 
contr61e 36, et N premiers modules additionneurs 40 A , 
5 40 N respectivement disposes sur les N voies de haut-parleur 
et recevant chacun en entr6e le signal sonore regu sur la 
voie consid6r&e et le signal de sortie du module de produc- 
tion de signaux auxiliaires associS S cette voie. 

L'annuleur d'6cho reprfesente sur la figure 4 

10 comprend en outre N x N modules de filtrage adaptatif 42 A j 
ou i et j sont des entiers variant de 1 a N, les N modules 

de filtrage adaptatif 42 A tl , , 42 i v 6tant associds d la 

i&ne voie de microphone, 1 < i £ N, et 6tant respectivement 
associSs aux N voies de haut-parleur, conune le montre la 

15 figure 4. Pour la ieme voie de microphone, 1 £ i £ N, les 
sorties des N modules de filtrage adaptatif 42 i#1 , 42 A N 

sont connect&es d un second module additionneur 44 1 , La 
sortie de chaque second module additionneur 44^ 1 £ i £ N, 
est connect6e & 1 1 entr6e soustractrice d 1 un module soustrac- 

20 teur 46 t recevant sur son entree non soustractrice le signal 
regu par le i&me microphone 32 x . La sortie de chaque module 
soustracteur 46 A commande 1' ensemble des N modules de 
filtrage adaptatif 42 t lf . 42 A N , par 1 ' intermediaire d'un 
module (non repr6sent6) permettant de calculer une erreur 

25 d' estimation & partir du resultat fourni par le module 
soustracteur 46^ et par 1 ' interm6diaire d 1 un module ( non 
repr6sent6) pour corriger de fagon it6rative les coeffi- 
cients des N modules de filtrage adaptatif 42 llf 42 1N , 
en fonction de 1' erreur d' estimation obtenue. 

30 La figure 5 reprfesente plus en detail la structure 

du sous-ensemble du module 34 de calcul de signaux auxiliai- 
res form6 par le module de contrOle 36 et les N modules 38 lf 
38 N de production de signaux auxiliaires. 

Le signal sonore regu sur une voie de haut-parleur 

35 donn6e est fourni en entr6e d un module 48 de formation de 
blocs. Le module 48 fournit en sortie des blocs contenant 
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chacun un nombre predetermine d 1 6chantillons du signal 
sonore regu. Deux blocs successifs se recouvrent mutuelle- 
ment sur un nombre predetermine d'6chantillons. Les exemples 
num6riques donnas dans la description d'un mode particulier 
5 de realisation du proc6d6 restent valables. 

Le module 48 de formation de blocs est connecte d un 
module de ponderation 50 qui applique a chaque bloc une 
fenStre d ' apodisation, par exemple du type fenetre de 
Hamming ou de Hanning. 

10 Le module de pond6ration 50 est connects & un module 

52 de calcul de transform6e de Fourier. A titre d' exemple 
non limitatif , le module 52 peut effectuer une transformee 
de Fourier rapide sur 512 points. 

Le module 34 de calcul de signaux auxiliaires 

15 comprend 6galement un module 54 de calcul des valeurs des 
frequences inf£rieures et superieures de bandes critiques 
adjacentes couvrant le domaine de frequences auquel appar- 
tient le signal sonore regu. En variante, le module 54 peut 
consister en une table de consultation contenant ces 

20 valeurs. 

Le module 54 est connecte & un module 56 de calcul 
de fonctions de repartition de l'energie. L f entr6e du module 
56 est egalement connect6e au module 52 de calcul de 
transformee de Fourier. Le module 56 calcule et additionne 

25 l'energie des diff6rentes raies spectrales appartenant & 
chaque bande critique en utilisant les valeurs des frequen- 
ces limites des bandes critiques fournies par le module 54. 
Le module 56 fournit en sortie, pour chaque bande critique, 
la valeur d'une fonction de repartition de l'6nergie. 

30 Une seconde table de consultation 58 contient des 

valeurs repr6sentant des fonctions d'fetalement basilaire 
associees chacune £ une bande critique donn6e. La definition 
de ces fonctions a 6t6 donn6e plus haut en relation avec un 
mode particulier de realisation du proc6d6. 

35 La seconde table de consultation 58 est connect6e & 

1* entr6e d r un module 60 permettant de calculer un spectre 
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6tal6 E t . L'entrfee du module 60 est 6galement connect6e au 
module 56 de calcul de fonctions de repartition de l'6ner- 
gie. Le module 60 calcule, pour une bande critique donnfee, 
le produit de convolution de la fonction de r6partition de 
5 l f 6nergie et de la fonction d'6talement, de fagon & fournir 
en sortie le spectre 6tal6 pour cette bande critique. 

Dans le mode particulier de realisation reprfesentfe 
sur la figure 5, le module 34 de calcul de signaux auxiliai- 
res comprend en outre un module 62 de calcul d' indices de 

10 tonalit6 a ± respectivement associ6s a des bandes critiques 
donn6es. Le module 62 est connects en sortie du module 52 de 
calcul de transformfee de Fourier et utilise les valeurs des 
6nergies des raies spectrales appartenant aux diff6rentes 
bandes critiques. 

15 Dans ce mode particulier de realisation, le module 

62 de calcul d f indices de tonalit6 est connect* S un module 
64 de calcul de facteurs de correction O t respectivement 
associ6s h des bandes critiques donnfees. Le module 64 
calcule les facteurs de correction Oi £ partir des indices 

20 de tonalite a A . 

Dans un mode particulier de realisation, le module 
62 de calcul d ? indices de tonalite fournit un indice de 
tonalite a t constant, identique pour toutes les bandes 
critiques, dfefini par : 

25 a 4 - a - min( SFM/SFM m , x , 1) 

ou min(a, b) d6signe la plus petite des valeurs a et b, 

ou SFM. M est un param6tre de valeur pr6d6termin6e en dB 

assocife 4 un signal sinusoidal pur, et 

ou SFM « 10 log 10 G/A, ou log 10 d6signe le logarithme decimal, 
30 G d6signe la moyenne g6om6trique de l'energie sur un nombre 
pr6d6termin6 de points de la transformee de Fourier et A 
d&signe sa moyenne arithmfetique sur le meme nombre de 
points. Dans ce mode de realisation, le module 64 de calcul 
de facteurs de correction 0i fournit en sortie pour chaque 
35 bande critique un facteur de correction 0 ± d6f ini par : 

Oi » max (SO, a L . (k x + B 4 ) + (1 - aj.k,) 
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oO max(a, b) designe la plus grande des valeurs a et b, 

ou SO, k x et k 2 sont des paramStres de valeur predetermine 

en dB, 

o£i a A est 1 f indice de tonality associ6 & la bande critique 
5 consid6r6e, et 

ou B L designe la frequence de la bande critique en Bark, 
A titre d'exemple non limitatif, les paramStres SFM.^, SO, 
k x et k 2 peuvent avoir les valeurs num6riques indiqu6es plus 
haut en relation avec un mode particulier de realisation du 
10 proc6d6. 

Le module 64 de calcul de facteurs de correction 0 L 
et le module 60 de calcul de spectres 6tal6s sont connect6s 
& l'entr6e d f un module 66 de calcul de seuils de masquage 
frequent iel T t respect ivement associ6s & des bandes criti- 
15 ques donnSes. Le module 66 fournit en sortie pour chaque 
bande critique un seuil de masquage fr6quentiel T £ d6fini 
comme suit : 

10 log 10 T x = 10 log 10 E t - 0 A - C t 
oil E t est le spectre 6tal6 et 0 4 le facteur de correction 

20 associ6s 4 cette bande critique, et ou C £ est un parametre 
suppl6mentaire de correction de seuil, 

Le paramStre C t est fourni par le module de contrdle 
36 mentionn6 plus haut, Les connexions du module 36 aux N 
voies de haut-parleur, representees sur la figure 4, ont 6t6 

25 omises sur la figure 5 pour plus de clart6. 

Dans un mode particulier de realisation, le module 
de contrdle 36 calcule le parametre supplemental re de 
correction de seuil C L pour chaque bande critique S partir 
du niveau respectif des signaux des N voies et d'un ensemble 

30 de valeurs de la fonction d f intercorr61ation d'au moins 
certaines paires de ces signaux, comme on l'a d6crit plus 
haut en relation avec un mode particulier de realisation du 
proc6d6. On a vu que dans ce cas, le paramfetre Cj est une 
fonction croissante des valeurs de la fonction d ' intercorr6- 

35 lation* 
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Dans un autre mode particulier de r6alisation, le 
module de contrSle 36 calcule le param^tre C L pour chaque 
bande critique & partir de la somme des Energies contenues 
dans les diverses bandes critiques de chaque voie et du 
rapport des Energies par bande critique d'au moins certaines 
paires des diffferentes voies, comme on l'a 6galement d6crit 
en relation avec un mode particulier de r6alisation du 
proc6d6. On a vu que dans ce cas, le paramfetre C t est 
d'autant plus important que la somme des Energies contenues 
dans les diverses bandes critiques est importante et que le 
rapport des energies par bande critique d'une paire de voies 
est peu different entre les diffferentes bandes critiques- 

Le module 34 de calcul de signaux auxiliaires 
comporte de plus un module 68 de g6n6ration de signaux 
synthfetiques. Le module produit, pour chaque voie de signal 
sonore, un signal ayant les caract6ristiques spectrales d'un 
bruit blanc et d6corr§l6 des signaux synth6tiques produit s 
pour les autres voies . 

La sortie du module 66 de calcul de seuils de 
masquage f rfequentiel est connectfee & un module de conversion 
fr6quentielle 69. Le module 69 effectue la conversion des 
seuils de masquage fournis par le module 66 de Bark en 
Hertz. 

Un module 70 regoit en entree les signaux synth6ti- 
ques produits par le module 68 et les seuils de masquage 
fr6quentiel calculfes par le module 66 et convertis par le 
module 69* Le module 70 am&ne le niveau de chaque signal 
synthfetique & la valeur du seuil de masquage frfequentiel 
associe, de fagon & obtenir un signal auxiliaire. 

La sortie du module 70 est connect6e a l'entrfee d'un 
module 72 de calcul de transformee de Fourier inverse. 

La sortie du module 72 est connect6e d l'entr6e d'un 
module de synthase 74. Le module de synthase 74 est connects 
& une mSmoire vive (non repr6sent6e) qui stocke le r6sultat 
du traitement des blocs d ! 6chantillons successifs. Le module 
de synthase 74 combine le rfesultat du traitement d'un bloc 
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avec le r6sultat du traitement du bloc pr6c6dent fourni par 
la m§ moire vive, par exemple en mettant en oeuvre un proc6d6 
de synthase par recouvrement et addition (en anglais 
" overlap- add " ) ou recouvrement et sauvegarde (en anglais 
5 "overlap -save" ) . 

Un module 76 de lecture s6quentielle est connects au 
module de synthase 74 et & la m6moire vive, et assure le 
passage de la cadence des donn&es de la cadence de blocs & 
une cadence d ' 6chanti lions* 

10 L' ensemble des modules 48 & 76 qui viennent d'fitre 

d6crits est equivalent 4 la s6rie de modules 38 x , . . . , 38 N de 
production de signaux auxiliaires represents de fagon 
simplififee sur la figure 4. En effet, chacun des modules 48 
d 76 contribue 4 produire un signal auxiliaire pour chaque 

15 voie de haut-parleur. 

Ainsi, la pr6sente invention permet d'injecter d 
chaque instant sur chaque voie de signal sonore regu un 
signal auxiliaire de niveau maximal, assurant efficacement 
1 ' accroissement de la vitesse de convergence, et dont le 

20 niveau est rendu par construction inf6rieur au seuil de 
masquage par le signal sonore re<?u. 

Le proc6d6 et le dispositif de la pr6sente invention 
s'appliquent 6galement en annulation d'echo mono-voie (cas 
N « 1), pour ameliorer la convergence du filtre adaptatif. 

25 lis sont applicables & tous les types d f annulation 

d'6cho acoustique, en particulier aussi bien dans le domaine 
temporel que dans le domaine fr&quentiel, ou encore en sous- 
bandes • 

Les figures 6a, 6b et 7a, 7b illustrent deux 
30 exemples de resultats comparatifs obtenus & l'aide du 
proc6d6 et du dispositif d f annulation d'6cho conformes & un 
mode particulier de realisation de la presente invention, 

Les graphes des figures 6a & 7b comprennent tous 
deux courbes. Les courbes en tirets concernent la voie 
35 stferfeophonique droite et les courbes en pointi!16s concer- 
nent la voie st6r6ophonique gauche. L'algorithme utilise 
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pour l f adaptation des coefficients des filtres adaptatifs 
est l'algorithme du gradient stochastique normalisfe ( NLMS ) • 
Les filtres adaptatifs ont une taille de 1024 coefficients. 
L'axe des abscisses repr6sente le nombre d'it6rations de 
5 l'algorithme d' adaptation, divis6 par un facteur 256. L'axe 
des ordonnees repr&sente l'6cart relatif, en dB, entre les 
r6ponses impulsionnelles r&elles des chemins de couplage 
acoustlque et les rgponses impulsionnelles des filtres 
adaptatifs d'un annuleur d'6cho st6r6ophonique. 

10 Le premier exemple comparatif, illustr6 par les 

figures 6a et 6b, correspond & un locuteur masculin. Le pas 
d' adaptation de l'algorithme NLMS est de 0,25. La figure 6a 
montre l'fecart relatif obtenu en utilisant la prdsente 
invention et la figure 6b montre l'6cart relatif obtenu en 

15 utilisant un proc6d€ et un dispositif classiques d T ambula- 
tion d'echo, ne faisant notamment pas intervenir de signaux 
auxiliaires. 

Le deuxidme exemple comparatif, illustr6 par les 
figures 7a et 7b, correspond a un locuteur f6minin. Le pas 

20 d' adaptation de l'algorithme NLMS est de 0,125. Les figures 
7a et 7b sont les analogues respectifs des figures 6a et 6b. 

Dans les deux exemples comparatif s, 1 ' utilisation de 
la prfesente invention assure toujours un 6cart relatif 
inferieur a 1 (c f est-i-dire negatif en dB), ce qui n'est pas 

25 le cas du procede et du dispositif classiques pour la voie 
droite (courbes en tirets). De plus, dans les cas d f utilisa- 
tion de la prdsente invention, l'ecart sur la voie gauche 
(courbes en pointill6s) est infferieur de plusieurs dB & 
l'6cart observe sur les courbes correspondant au proc6d6 et 

30 au dispositif classiques. La presente invention prfesente 
done des performances supferieures, permettant de garantir un 
bon comportement de poursuite de l f annuleur d'6cho en 
presence de variations des chemins de couplage acoustlque. 

Par ailleurs, dans les deux exemples, la correction 

35 appliqu&e aux seuils de masquage est de 10 dB. On pourrait 
tolferer une r6duction de quelques dB de cette correction, ce 
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qui amfeliorerait encore les performances du proc6d6 et du 
dispositif propos6s, alors que la parole pergue, bien que 
16g§rement modifi6e par le bruit ajout6, conserverait une 
qualitfe tout & fait acceptable dans une application du type 
teleconference . 
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REVEND I CATIONS 

1. Procdd6 d'annulation d'6cho sur N voles de signal 
sonore ayant chacune un haut-parleur et un microphone 
5 associ6s, N 6tant un entier sup&rieur ou 6gal & 1, suivant 
lequel, sur chacune des N voies : 

(a) on engendre un signal synth6tique ayant les 
caract6ristiques spectrales d'un bruit blanc, le spectre de 
ce signal s'&tendant sur plusieurs bandes de frequence 

10 adjacentes, et ce signal synth6tique 6tant d6corr616 des 
signaux synth6tiques engendr6s sur les autres voies ; 

(b) pour chaque bande de frequence, on calcule un 
seuil de masquage frSquentiel correspondant au signal 
associg au haut-parleur de la voie consid6r6e en utilisant 

15 des propri6t6s de la perception auditive humaine ; 

(c) dans chaque bande de frequence, on amdne le 
niveau du signal synthfetique d la valeur du seuil de 
masquage frfequent iel assocife, de fagon d obtenir un signal 
auxiliaire ; 

20 (d) on additionne le signal auxiliaire au signal 

associe au haut-parleur de la voie consider6e, les seuils de 
masquage frequentiel ayant 6t6 precedemment calcul6s de 
fagon a rendre le signal auxiliaire inaudible, et les 
signaux auxiliaires des N voies fetant mutuellement d6corr6- 

25 16s ; 

(e) on fournit le signal obtenu pr6c6demment en 
entr6e & un filtre adaptatif dont les coefficients forment 
une estimation de la rfeponse impulsionnelle du chemin de 
couplage acoustique entre le haut-parleur et le microphone 

30 associ6s & la voie de signal sonore consid6r6e ; 

( f ) on additionne les signaux obtenus en sortie des 
N filtres adaptatifs associSs & chaque voie de microphone et 
on soustrait le signal resultant du signal regu par le 
microphone associe & cette voie ; 
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(g) on calcule une erreur d ' estimation 4 partir de 
la difference obtenue £ 1' issue de la soustraction pr6c6- 
dente ; 

(h) on corrige les coefficients des f litres adapta- 
5 tifs associ6s ft la voie de microphone consid6r6e en fonction 

de 1 ' erreur d' estimation associ6e. 

2. Proc6d6 selon la revendication 1, caract6ris6 en 
ce que, pour calculer chaque seuil de masquage fr6quentiel, 
(bl) on forme des blocs contenant chacun un nombre 
10 predetermine d f 6chanti lions de signal sonore, deux blocs 
success! fs se recouvrant mutuellement sur un nombre prede- 
termine d ' 6chanti lions ; 

(b2) on pondere les echantillons de chaque bloc par 
une fenetre d' apodisation ; 
15 (b3) on calcule la transform6e de Fourier de chaque 

bloc ; 

(b4) on divise le domaine des frequences en plu- 
sieurs bandes critiques adjacentes ayant des frequences 
inf6rieures et sup6rieures determin6es ; 
20 puis, pour chaque bloc et dans chaque bande critique : 

(b5) on additionne I'energie des diff6rentes raies 
spectrales appartenant a la bande critique, de fagon & 
obtenir la valeur, pour cette bande critique, d'une fonction 
de repartition de I'energie ; 
25 (b6) on calcule le produit de convolution de la 

fonction de repartition de l'6nergie et d'une fonction 
d'6talement basilaire, obtenue & partir d'une table de 
consultation, de fagon S obtenir un spectre 6tal6 E ± ; 

(b7) on calcule un indice de tonalite a L d partir des 
30 energies des raies spectrales appartenant aux differentes 
bandes critiques ; 

(b8) on calcule un facteur de correction 0 t k partir 
de 1' indice de tonalite calcule pr6c6demment ; 

(b9) on calcule un seuil de masquage fr6quentiel T A 
35 a partir du spectre 6tal6 E i et du facteur de correction 0 4 , 
d6fini comme suit : 
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10 log 10 T t = 10 log 10 E t - 0 A - C t 
ou C t est un paramfetre supplfementaire de correction de seuil 
calculfe k partir de plusieurs caracteristiques des signaux 
sonores des N voles* 
5 3. Procfedfe selon la revendication 1 ou 2, caractfe- 

risfe en ce que les caractferistiques servant & calculer ledlt 
paramfetre supplfementaire de correction de seuil C t compren- 
nent le niveau respect if des signaux des N voles et un 
ensemble de valeurs de la fonction d ' intercorrfelation d'au 

10 moins certaines paires de ces signaux, ledit paramfetre 
supplfementaire de correction C t fetant une fonction crois- 
sante des valeurs de cette fonction d' intercorrfelation. 

4. Procfedfe selon la revendication 1 ou 2, caractfe- 
risfe en ce que les caractferistiques servant 4 calculer ledit 

15 paramfetre supplfementaire de correction de seuil C A compren- 
nent la somme des Energies contenues dans les diverses 
bandes critiques de chaque voie et le rapport des Energies 
par bande critique d'au moins certaines paires des differen- 
tes voies, ledit paramfetre supplfementaire de correction C L 

20 fetant d ' autant plus important que la somme des fenergies 
contenues dans les diverses bandes critiques est importante 
et que le rapport des fenergies par bande critique d'une 
paire de voies est peu different entre les diffferentes 
bandes critiques* 

25 5. Annuleur d'echo £ N voies de signal sonore ayant 

chacune un haut-parleur (30 : , . .., 30 N ) et un microphone 
(32 x/ 32 N ) associfes, N fetant un entier supferieur ou fegal 

d 1. comprenant : 

- des moyens (34) de calcul de signaux auxiliaires, 

30 comportant : 

des moyens ( 48 ) de formation de blocs conte- 
nant chacun un nombre prfedfeterminfe d ' fechantillons du 
signal sonore, deux blocs successifs se recouvrant 
mutuellement sur un nombre prfedfeterminfe d'fechantil- 
35 Ions ; 
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des moyens (50) de ponderation des 6chantil- 
lons de chaque bloc par une fenetre d f apodisation, 
disposes en sortie des moyens de formation de 
blocs ; 

5 des moyens (52) de calcul de transform6e de 

Fourier, disposes en sortie des moyens de pondera- 
tion ; 

des moyens (54) de calcul des valeurs de 
frequences inf6rieures et superieures de bandes 
10 critiques adjacentes, ou une table de consultation 

contenant ces valeurs ; 

des moyens (56) pour calculer une fonction de 
repartition de l'6nergie sur les bandes critiques, 
la valeur, dans une bande critique donn6e, de la 
15 fonction de repartition de l ? 6nergie 6tant calcul6e 

en additionnant l'energie des diff6rentes raies 
spectrales appartenant & cette bande critique ; 

une seconde table de consultation (58), 
contenant des valeurs repr6sentant des fonctions 
20 d'6talement basilaire associ6es chacune A une bande 

critique donn6e ; 

des moyens ( 60 ) pour calculer un spectre 6tal6 
E L , en calculant, pour une bande critique donn6e, le 
produit de convolution de la fonction de 
25 repartition de 1'energie et de la fonction d'6tale- 

ment associfees & cette bande critique ; 

des moyens (64) de calcul d'un facteur de 
correction 0i associe & une bande critique donn6e ; 

des moyens (66) de calcul d'un seuil de 
30 masquage fr6quentiel T t associe & une bande 

critique donn6e, 3t partir du spectre 6tal6 et du 
facteur de correction 0 t associes a cette bande 
critique, le seuil de masquage frequentiel T L etant 
d6fini comme suit : 
35 10 log 10 T t - 10 log 10 E t - 0 t - C t 
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ou C x est un param6tre supplementaire de correction 
de seuil ; 

des moyens (68) de g6n6ration de signaux 
synthetiques ayant les caract6ristiques spectrales 
5 d'un bruit blanc ; 

des moyens (36) de controle, connect6s aux N 
voies de signal sonore, pour calculer pour chaque 
bande critique ledit param6tre suppl6mentaire de 
correction de seuil Ci & partir de plusieurs carac- 
10 t6ristiques des signaux sonores des N voies ; 

des moyens ( 70 ) pour amener le niveau d ' un 
signal synthfetique & la valeur du seuil de masquage 
frequent iel associ6, de fagon a obtenir un signal 
auxiliaire ; 

15 des moyens (72) de calcul de transform6e de 

Fourier inverse ; 

des moyens de m6moire pour stocker le rdsultat 
du traitement des blocs successifs ; 

des moyens (74) pour combiner le r6sultat du 
20 traitement d'un bloc avec le r6sultat du traitement 

du bloc pr6c6dent ; 

des moyens (76) de lecture s6quentielle relifes 
auxdits moyens de m6moire, pour faire passer la 
cadence des donnfees de la cadence bloc S la cadence 
25 6chantillon ; 

N premiers moyens additionneurs (40 X/ 
40 H ) respectivement dispos6s sur les N voies de 
signal sonore, pour additionner sur chaque voie le 
signal sonore regu et le signal auxiliaire associ£ ; 
30 - N x N moyens de filtrage adaptatif (42 t lf . • . , 42 L M , 

. 42 M lf . 42 NN ), comprenant N moyens de filtrage 

adaptatif associ6s S chacune des N voies de microphone, et 
ayant chacun des coefficients qui forment une estimation de 
la r6ponse impulsionnelle du chemin de couplage acoustique 
35 entre le haut-parleur et le microphone associ£s & une des N 
voies de signal sonore ; 
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- des seconds moyens additionneurs (44 1# 44 N ), 
connectes en sortie des N moyens de filtrage adaptatif 
associes & chaque voie de microphone, pour additionner les 
signaux de sortie de ces N moyens de filtrage adaptatif ; 

5 - N moyens soustracteurs (46^ . 46 M ), disposes 

chacun en sortie desdits seconds moyens additionneurs (44 w 
44 N ) et connectes respect ivement aux microphones des N 
voies de signal sonore, pour soustraire sur chaque voie de 
microphone le signal obtenu en sortie desdits seconds moyens 
10 additionneurs (44 x , 44 N ) du signal regu par le micro- 

phone de cette voie ; 

- des moyens pour calculer sur chaque voie de 
microphone une erreur d f estimation & partir du r6sult:at 
fourni par le moyen soustracteur connects au microphone de 

15 cette voie ; 

- des moyens pour corriger de fagon it6rative les 
coefficients respectifs des N moyens de filtrage adaptatif 
associ&s d chacune des N voies de microphone, en fonction de 
1 ' erreur d' estimation associ&e & chaque voie de microphone. 

20 6. Dispositif selon la revendication 5, caract&ris6 

en ce que les caracteristiques servant S calculer ledit 
paramfetre suppl6mentaire de correction de seuil C A compren- 
nent le niveau respectif des signaux des N voies et un 
ensemble de valeurs de la fonction d 1 intercorrfelation d'au 

25 moins certaines paires de ces signaux, ledit paramStre 
supplfementaire de correction C t &tant une fonction crois- 
sante des valeurs de cette fonction d 1 intercorr61ation. 

7. Dispositif selon la revendication 5, caract6ris6 
en ce que les caract6ristiques servant & calculer ledit 

30 paramfetre suppl6mentaire de correction de seuil C A compren- 
nent la somme des Energies contenues dans les diverses 
bandes critiques de chaque voie et le rapport des Energies 
par bande critique d'au moins certaines paires des diff6ren- 
tes voies, ledit paramdtre suppl6mentaire de correction C x 

35 §tant d*autant plus important que la somme des Energies 
contenues dans les diverses bandes critiques est importante 
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et que le rapport des Energies par bande critique d'une 
paire de voies est peu different entre les diff6rentes 
bandes critiques. 

8. Dispositif selon la revendication 5, caract6ris6 
5 en ce qu'il comprend en outre des moyens (62) de calcul d'un 

indice de tonalite a L a partir des Energies des raies 
spectrales appartenant aux diff6rentes bandes critiques, 
lesdits moyens ( 64 ) de calcul du f acteur de correction 0 t 
calculant ledit facteur de correction 0 1 & partir dudit 
10 indice de tonalit6 a t . 

9. Dispositif selon la revendication 8, caract6ris6 
en ce que lesdits moyens (64) de calcul du facteur de 
correction 0i fournissent un facteur de correction 0 t d6fini 
par : 

15 0 t - max (SO, (kj + B L ) + (1 - a t ) .k 2 ) 

ou max(a, b) d6signe la plus grande des valeurs a et b, 

o£i SO, k x et k 2 sont des param&tres de valeur pr6d6termin6e 

en dB, 

oCi a A est 1 ' indice de tonalit6 associe a la bande critique 
20 considferee, et 

ou B L d6signe la frequence de la bande critique en Bark* 

10. Dispositif selon la revendication 9, caract6ris6 
en ce que SO = 24,5 dB, k x = 14,5 dB et k 2 * 5,5 dB. 

11. Dispositif selon la revendication 8, 9 ou 10, 
25 caract6ris6 en ce que les moyens (62) de calcul de 1 ' indice 

de tonalit6 fournissent un indice de tonalit§ constant, 
identique pour toutes les bandes critiques et defini par : 

a t = a = min(SFM/SFM >l „, 1) 
oil min(a, b) d&signe la plus petite des valeurs a et b, 
30 ou SFM^ est un paramdtre de valeur pr6d6termin6e en dB 
assocife a un signal sinusoidal pur, et 

ou SFM = 10 log 10 G/A, ou log 10 d6signe le logarithme d6cimal, 
G dfesigne la moyenne geomStrique de l^nergie sur un nombre 
pr6d6termin6 de points de la transform6e de Fourier et A 
35 d6signe sa moyenne arithm6tique sur le meme nombre de 
points. 
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12. Dispositif selon la revendication 11, caracte- 
rise en ce que SFM B „ = - 60 dB. 
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